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The Rashba effect is one of the most striking manifestations of spin-orbit coupling in solids, 
and provides a cornerstone for the burgeoning field of semiconductor spintronics. It is typically 
assumed to manifest as a momentum-dependent splitting of a single initially spin-degenerate band 
into two branches with opposite spin polarisation. Here, combining polarisation-dependent and 
resonant angle-resolved photoemission measurements with density-functional theory calculations, 
we show that the two “spin-split” branches of the model giant Rashba system BiTel additionally 
develop disparate orbital textures, each of which is coupled to a distinct spin configuration. This 
necessitates a re-interpretation of spin splitting in Rashba-like systems, and opens new possibilities 
for controlling spin polarisation through the orbital sector. 


INTRODUCTION 

The ability to controllably engineer spin splittings of 
electronic states is a key goal in the search for spin- 
tronic materials [lj. A particularly successful strategy 
has been the lifting of spin degeneracy via spin-orbit 
coupling in the presence of structural inversion asym¬ 
metry. Termed the Rashba or Bychov-Rashba effect [2] 
this manifests through a spin-momentum locking of the 
quasiparticles, stabilising a pair of Fermi surfaces which 
are typically assumed to exhibit counter-rotating chiral 
spin textures [3 r The ability to electrostatically con¬ 
trol the strength of such spin splitting 0H7] has led to 
prospects for all-electrical manipulation of electron spin 
precession [8], offering new prototypical schemes of semi¬ 
conductor spintronics p] . The quest to practically realise 
such devices, and to operate them without cryogenic cool¬ 
ing, has motivated a major search for materials which 
can host stronger spin splittings than can typically be 
achieved in conventional semiconductors [6j I9HT2] . 

A giant Rashba-like spin splitting was recently dis¬ 
covered for bulk conduction and valence band states of 
the bismuth tellurohalide semiconductors H3HH3. Aris¬ 
ing due to a combination of bulk inversion asymmetry 
(Fig. Em strong atomic spin-orbit coupling, and a neg¬ 
ative crystal field splitting of the valence bands na, 
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Rashba parameters amongst the highest of any known 
material have been uncovered, together with a counter¬ 
rotating chiral Fermi surface spin texture na esj. 
Exploiting element- and orbital-selective angle-resolved 
photoemission (ARPES), we show here that a complex 
interplay between atomic, orbital, and spin degrees of 
freedom significantly enriches this picture. We expect 
our findings to be broadly applicable across other strong 
spin-orbit Rashba systems. 


RESULTS 

First, we summarise the generic electronic structure of 
the Te-terminated surface of BiTel (Fig. [l]) . Previous 
measurements have shown this termination to support a 
near-surface electron accumulation nim. We clearly 
observe two 2D subbands (SB1 and SB2 in Fig.[l|c), see 
also Supplementary Fig. SI) formed within the result¬ 
ing quantum well. Each subband hosts two branches 
with a separation that grows approximately linearly 
with momentum away from their crossing at k pa 0. 
This is a hallmark of Rashba-like spin splitting. We 
extract a large Rashba energy Er = 120 db 10 meV 
(85 d= 5 meV) and momentum offset of the band bot- 
tom k 0 = 0.055 ± 0.005 A " 1 (0.050 ± 0.005 A" 1 ) for the 
first (second) subband, respectively, supporting previous 
studies which established this material as a model host 
of giant spin splittings mmm- 
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FIG. 1. Surface electronic structure of BiTel. (a) 

ARPES measurements of the Fermi surface and near-E f band 
dispersions measured (hv — 52 eV, p-polarisation) from the 
Te-terminated surface of BiTel. A lack of inversion symmetry 
of the bulk crystal structure (b) together with strong spin- 
orbit coupling mediates a large Rashba-like spin splitting. 
Additionally, the polar nature of the Te-terminated surface 
induces a strong downward band bending, causing a ladder 
of Rashba-split subband (SB) states to emerge in the near¬ 
surface quantum well, (c) These are clearly resolved in mea¬ 
surements of the dispersion along T — M. The conventional 
spin texture associated with such Rashba splitting is shown 
schematically by coloured arrows, with the spin expected to 
lie predominantly in the surface plane. 


While typically treated in a single-band picture, the 
electronic wavefunction for each branch of a Rashba- 
split state can more generally be written as = 
r a c i r'tPhr, where, following the notation in Ref. [18!, 
i is the atomic index, r G {p x ,Py,Pz} and a are the or¬ 
bital and spin indices, respectively, ^ >r are atomic wave- 
functions, and <? iT are complex coefficients. Neglecting 
spin-orbit coupling, our calculations predict a conduction 
band in BiTel predominantly derived from Bi p z orbitals 
(see Supplementary Fig. S2). Including such effects, how¬ 
ever, not only permits it to become strongly spin-split via 
Rashba-like interactions, but also promotes significant 
orbital mixing. In general, therefore, multiple c? r can 
be expected to become non-negligible. For a complete 
description of the Rashba-split states, it is therefore es¬ 
sential to consider the interplay of the underlying atomic, 


orbital, and spin components. To disentangle these, we 
combine two powerful features of ARPES: characteristic 
selection rules for photoemission using linearly-polarised 
light, allowing us to directly probe the orbital wavefunc¬ 
tion B21ID], and resonant photoemission to provide ele¬ 
mental sensitivity [21]. Such resonant enhancements are 
evident in Fig. [2j[a,b) . They cause the spectral weight of 
the conduction band states to become strongly peaked 
at photon energies around 26 and 28 eV, close in energy 
to the binding energy of the Bi 5d 5 / 2j 3/2 core levels, with 
functional forms that are well described by Fano line- 
shapes. This points to a significant Bi-derived atomic 
character of the lowest conduction band states, consistent 
with theoretical calculations [22[. We exploit this, selec¬ 
tively probing “on-resonance” to unveil the Bi-projected 
spectral function. 

The resulting measurements of the dispersion reveal 
pronounced momentum-dependent spectral weight vari¬ 
ations (Fig. j2jc)). We focus on the first subband 
(SB1), which is most clearly visible across our measure¬ 
ments. When probed using p-polarised light, we find 
stronger spectral weight for the inner branch of this sub¬ 
band, while the outer branch is significantly more pro¬ 
nounced when measured using 5-polarisation. Selection 
rules [19] dictate that, of the p-orbitals, the former mea¬ 
surement should be sensitive to p z - and p^-derived or¬ 
bital character, while the transition matrix element is 
only non-vanishing for photoemission from p y orbitals in 
the latter case (our measurement geometry is shown in 
Fig. [3][a)). The asymmetric spectral weight distributions 
within and between these measurements immediately es¬ 
tablishes that the two spin-split branches of the disper¬ 
sion host a markedly different orbital makeup. Moreover, 
when measuring with a photon energy only 2 eV higher 
(Fig- id)), we find an almost complete reversal of these 
matrix element variations, with greater spectral weight 
for the inner branch of the lowest subband when measur¬ 
ing using 5-polarised light. No longer on resonance, these 
measurements will not a priori be dominated by the Bi- 
derived spectral weight, and we therefore conclude that 
the orbital textures of the Rashba-split states are also 
strongly dependent on their underlying atomic charac¬ 
ter, as discussed in detail below. 

Fig. §b, c) shows constant binding energy surfaces 
(CESs) measured on-resonance with 5-polarised light, to 
selectively probe the momentum-space distribution of the 
Bi p y orbital character. The outer band of each CES ex¬ 
hibits strong spectral weight along our measurement di¬ 
rection (at positive and negative k y ), with a pronounced 
suppression at positive and negative k x , over the entire 
occupied bandwidth of the conduction bands. This re¬ 
mains qualitatively unchanged under rotation of the az¬ 
imuthal orientation of the sample. Such behaviour is 
indicative of a p y -like orbital aligned along our measure¬ 
ment direction, irrespective of whether the T—M or T—K 
crystallographic direction is oriented to this jl8j [20;. The 
in-plane Bi-projected charge density must therefore be 
oriented such that its lobes point predominantly radially 
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FIG. 2. Disentangling intertwined atomic and orbital characters, (a) Fermi level momentum distribution curve 
(. Ef =L 15 meV) measured along V — M as a function of probing photon energy using p-polarised light. No dispersion is 
observed, indicating two-dimensional states consistent with our assignment of quantum well subbands, while strong matrix 
element variations give rise to pronounced intensity modulations. These are shown for the outermost band (kp ~ —0.2 A -1 ) 
as a function of binding energy in (b), revealing characteristic intensity enhancement due to resonant photoemission at the 
Bi O-edge. (c,d) Corresponding dispersions (along Y — M ) measured on- (his — 28 eV) and off- (his = 30 eV) resonance, 
respectively, with p-polarised (left) and s-polarised (centre) light. The difference between these (right) indicates a band and 
elemental-dependent orbital polarisation. 


outwards from the outer CES. 

A similar conclusion can be drawn for the inner band 
of the CES at energies below the Dirac point formed by 
the crossing of the two spin-split branches of the disper¬ 
sion. As already evident 100 meV below the Dirac point 
in Fig.[3jb, c), and more clearly shown in the angular de¬ 
pendence of spectral weight extracted around the CES in 
Fig. B a), the spectral weight is still peaked along the k y 
direction (a = 0 and it in Fig. |4ja)), with a suppression 
of spectral weight along k x (a = it/2). Indeed, while 
the band edge turning point in the dispersion leads to 
a van Hove singularity in the density of states [23] the 
inner band of the CES below the Dirac point smoothly 
evolves into the outer band as it moves through this turn¬ 
ing point, and thus hosts qualitatively the same orbital 
texture. 

Moving up through the Dirac point, however, the angu¬ 
lar distribution of spectral weight around the inner CES 
gradually flattens and then inverts to become peaked 
at a = 7r/2. This points to a tangential alignment of 
the in-plane Bi-derived orbital character above the Dirac 
point, with stronger p y character at positive and nega¬ 
tive k x than positive and negative k y . This behaviour is 
well-reproduced by our first-principles calculations. As 


shown in Fig. [4[c), the in-plane p y : p x orbital polari¬ 
sation around the calculated CESs of the lower subband 
show a clear peak at a = n/2 above the Dirac point, 
and 0 and i r below. To capture this generic behaviour, 
we extract a relative spectral weight angular distribution 
factor from our measurements : 

a / \ _ C* = 7r/2(^) — A* = 0,7T (^0 / -| \ 

where I a ( uj) is the spectral weight at angle a and energy 
oj. This provides an experimental measure of the calcu¬ 
lated “orbital polarisation ratio” used to describe topo¬ 
logical states in Ref. m- Here, this parameterises the 
relative strength of spectral weight, and thus p y orbital 
character, along the k x and k y directions of our mea¬ 
sured CESs, reflecting the strength (|A|) and alignment 
of radial (A < 0) and tangential (A > 0) in-plane orbital 
textures. As shown in Fig. [4^b), we find this is positive 
at energies above the Dirac point for the on-resonance 
measurements, consistent with the dominantly tangen¬ 
tial in-plane Bi-projected charge density distribution of 
the inner band assigned above. With decreasing energy 
towards the Dirac point, A is suppressed to zero, indi¬ 
cating a loss of orbital polarisation, before growing again 
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FIG. 3. ^-projected spectral weight distribution, (a) 

Experimental geometry for our measurements, and (b-d) re¬ 
sulting constant binding energy surfaces measured using s- 
polarised light to probe the p y orbital character at 200 meV 
above, exactly at, and at 100 meV below the Dirac point (Ed) 
formed by the crossing of the two spin-split branches of the 
lowest subband. On-resonance measurements (hv — 28 eV) 
with the scattering plane aligned to (b) T — M and (c) 
r — K, and (d) off-resonance ( his = 30 eV, scattering plane 
along T — K) measurements show pronounced angular varia¬ 
tions in spectral weight, characteristic of strongly momentum- 
dependent orbital textures. 


but with opposite sign below the Dirac point. Thus, we 
have experimentally observed a gradual crossover from a 
radial to a tangential alignment of the in-plane Bi orbital 
character with the switch occurring exactly at the Dirac 
point within our experimental error. This is strikingly 
similar to a switch in orbital texture that was recently 
reported at the Dirac point of the spin-polarised surface 
states of topological insulators [T8l 120] . Our observation 
of such a crossover in a topologically-trivial compound 
establishes this as a general feature of strongly spin-orbit 
coupled systems, a point we return to below. 

Intriguingly, extracting A (a;) from our measurements 
performed off-resonance, we find an opposite trend 
(Fig- gb)). Above the Dirac point, A < 0 indicating a 
radial alignment of the in-plane orbital-projected charge 
density, as also clearly evident in the spectral weight dis¬ 
tribution visible in Fig. §d). This is again smoothly 
suppressed to zero approaching the Dirac point before be¬ 


coming positive, revealing a tangential orbital configura¬ 
tion, at energies below Dirac point. We attribute this as 
a signature of the in-plane Te orbital polarisation, which 
well matches that found in our calculations (Fig. Qd)). 
It can be seen in our measurements here at 30 eV photon 
energy as, with no Bi resonant enhancement, the photoe¬ 
mission cross-section is higher for Te 5 p than Bi 6p or¬ 
bitals m and, with Te being located right at the surface, 
there is no depth-dependent attenuation due to the sur¬ 
face sensitivity of photoemission. Our calculations (Sup¬ 
plementary Fig. S2) reveal a strong in-plane Te char¬ 
acter of the inner branch of the dispersion, which can 
thus dominate the spectral weight of the inner branch of 
the CES in these off-resonance measurements. For the 
outer branch, however, there is little calculated in-plane 
Te weight (Supplemental Fig. S2), and so the weak spec¬ 
tral features visible for the outer branch in Fig. J^d) is 
still reflective of the Bi in-plane orbital texture. 


DISCUSSION 


Together, these measurements and calculations reveal 
that spin-orbit coupling induces a complex atomic- and 
momentum-dependent hierarchy of orbitally-polarised 
components of the underlying electronic structure in 
BiTel, summarised schematically in Fig. [5] Our calcu¬ 
lations additionally reveal how each orbital component 
is, in turn, coupled to a disparate spin texture. We il¬ 
lustrate this for the Bi-derived states in Fig. [5| b-d); the 
Te-derived component is additionally shown in Supple¬ 
mentary Fig. S2. The in-plane spin-texture (S X:V ) pro¬ 
jected onto Bi p z orbitals yields a conventional counter¬ 
rotating chiral spin texture of neighbouring CESs at en¬ 
ergies above the Dirac point, characteristic of classic 
Rashba systems 0 0 and indeed observed experimen¬ 
tally for BiTel 13. 16]. In contrast, the spin texture 
is significantly more complex for the p x and p y orbital 
projection, with the in-plane spin component switching 
between tangential and radial around the CES. This re¬ 
sults from a coupling of the spin to the characteristic 
orbital textures, and is similar to that found recently for 
topological surface states |25H27] . 

For example, combined with our uncovered tangential 
(radial) in-plane Bi orbital texture above (below) the 
DP, this leads to a net clockwise spin rotation for the 
in-plane orbital projection of both the inner and outer 
CESs. This is in strong departure from the conven¬ 
tional picture of Rashba-like spin polarisation. Such in¬ 
plane orbital textures likely also underlie the unconven¬ 
tional spin topology predicted for certain p x , y -derived 
states in Bi/Cu(lll) [28] and Pb/Cu(lll) [2S] surface 
alloys. They are also broadly consistent with recent 
first-principles calculations that suggest the Rashba pa¬ 
rameter can be orbital-dependent in bismuth telluro- 
halides [30] . although here we reveal how orbital mixing 
can lock the spin components on the nominally “spin- 
split” CESs together, stabilising pronounced components 
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FIG. 4. Mapping the angle-dependent orbital wavefunctions. (a) Angular distribution of the Bi ^-projected spectral 
weight distribution measured on-resonance around the inner band of CESs (see e.g. Fig. K3[b)), each normalised to its average 
value. Above (below) the Dirac point, Ed , this is peaked at an azimuthal angle of a = tt/2 (0 and i r) indicative of a tangential 
(radial) in-plane orbital alignment, (b) The smooth evolution between these two configurations is captured by a relative 
spectral weight angular distribution factor (ADF) A (a;) — [I a = tt/ 2 (^) — 7 a= o,7r(^)]/[7 a=7r /2(k;) + Ia=o,-rr(w)\ which crosses zero 
at the Dirac point within our experimental error, (c) This behaviour is fully captured by the momentum-dependent Bi p y : p x 
in-plane orbital polarisation, £(u;,k), extracted from our density-functional calculations, (d) The in-plane orbital polarisation 
is reversed for the Te-projected component, with radial alignment above the Dirac point and tangential alignment below Ed , 
as also captured by the opposite ADF measured off-resonance (hv — 30 eV, see Fig. [3^c) ), summarised in (b). 


of the underlying wavefunction that host markedly non- 
Rashba-like spin textures in this system. Similar consid¬ 
erations hold for the Te-derived components, although 
with additional variations in the magnitude of the spin 
components projected onto out-of-plane vs. in-plane or¬ 
bital components between the CESs due to a greater out 
vs. in-plane orbital polarisation for Te (Supplementary 
Figs. S2 and S3). With each orbital component locked 
to a different spin texture, the fundamental requirement 
from time-reversal symmetry of spin-degeneracy at the 
Kramer’s point (the Dirac point formed in this system 
at k = 0) ensures equal contribution of in-plane p x and 
p y orbitals at this point. Thus, the vanishing of orbital 
polarisation at the Dirac point observed experimentally 
above can be simply viewed as an orbital analogue of 
Kramer’s spin degeneracy in time-reversal symmetric sys¬ 
tems. 

Away from the Dirac point, our calculations addition¬ 
ally predict a strong canting of the spin out of the sur¬ 
face plane for the Bi-derived orbitals. This grows in 
magnitude with increasing energy away from the Dirac 
point, where the CESs become increasingly hexagonally 
warped m We show evidence for this through circular 
dichroism measurements. We attribute such dichroism 
as a signature of unquenched orbital angular momen¬ 
tum m, which our calculations reveal is large and locked 
approximately opposite to the spin angular momentum 
due to the strong spin-orbit coupling. Circular dichroism 
is known to show a complex dependence on photon en¬ 


ergy in this system [32] as well as in other layered com¬ 
pounds such as Bi 2 Te 3 (Ref. [33]) . which can naturally 
arise as a consequence of inter-layer photoelectron inter¬ 
ference [H E2- To simplify such effects, we again 
perform our measurements on-resonance, selectively en¬ 
hancing states of Bi character (Fig. |5jf)). We find that 
the outer (hexagonally-warped) band in such measure¬ 
ments develops a pronounced six-fold modulation, which 
has previously been shown to reflect out-of-plane spin- 
canting in topological insulators [7, 34, .35 . Here, this 
provides the first experimental evidence for deviations 
from simple in-plane chiral spin textures in BiTel. 

As evident in Fig. |5jb-d) , (S z ) is substantially larger 
when projected onto the in-plane orbital components 
than for the p z projection. The emergence of a large out- 
of-plane spin canting in this system is thus intricately tied 
to the development of the in-plane orbital texture away 
from the Dirac point observed here. Together, our mea¬ 
surements and calculations establish a powerful role of 
in-plane atomic orbitals shaping the spin structure of the 
model Rashba compound BiTel, revealing a complex in¬ 
terplay between atomic makeup, anisotropic orbital tex¬ 
tures, and spin-momentum locking. While small quanti¬ 
tative variations may occur due to near-surface potential 
contributions to the Rashba spin splitting and possible 
surface-induced orbital reconstructions here, our findings 
should be broadly applicable to the spin-split bulk elec¬ 
tronic structure of BiTel as well as being generic to other 
strong spin-orbit Rashba systems, suggesting new routes 
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FIG. 5. Hierarchy of atomic, spin, and orbital coupling in BiTel. (a) Schematic representation of disparate in-plane 
orbital textures uncovered for the inner and outer “spin-split” branches a model Rashba system. Corresponding spin textures 
calculated from density-functional theory for Bi (b) p x , (c) p y , and (d) p z projections of CESs 200 meV above the Dirac point. 
The in-plane spin texture is shown by arrows and the out-of-plane by the background colour, both in units of /i/2. (e) Schematic 
of the in-plane spin textures coupled to the net in-plane and out-of-plane orbital textures, (f) Circular dichroism measurements 
performed on-resonance (Jw = 28 eV, CES at Ed + 200 meV), revealing a significant six-fold modulation for the outer band 
indicative of pronounced out-of-plane spin canting. 


to control spin splitting through the orbital sector. For 
example, exploiting structure-property relations to tune 
the competition of atomic spin-orbit coupling and crys¬ 
tal field splitting will allow control over the ratio of in¬ 
plane and out-of-plane orbital polarisation, thus mod¬ 
ulating the degree of non-Rashba-like spin components 
and out-of-plane spin canting. Together, this promises 
new prospects for the targeted design of optimised spin- 
tronic materials. 


MATERIALS AND METHODS 

Angle-resolved photoemission: ARPES measure¬ 
ments were performed at the CASSIOPEE beamline of 
SOLEIL synchrotron (France) and 105 beamline of Dia¬ 
mond Light Source (UK). Single-crystal samples of BiTel, 
grown by chemical vapour transport, were cleaved in-situ 
at the measurement temperature of 10 K. Measurements 
were performed using linear-horizontal, linear-vertical, 
and circularly-polarised synchrotron light at the photon 
energies described in the text, and reproduced on multi¬ 
ple samples. Scienta R4000 hemispherical electron anal¬ 
ysers were employed, with a vertical entrance slit and 
the light incident in the horizontal plane as shown in 
Fig-! a). BiTel is known to have domains of mixed Te 
and I termination, supporting surface electron and hole 
accumulation layers, respectively m- By monitoring the 
relative spectral weight of electron- and hole-like bands 
crossing the Fermi level and from characteristic core-level 


shifts in x-ray photoemission spectra, we aligned the syn¬ 
chrotron light spot on Te-terminated domains of the sam¬ 
ple. From core-level spectra, we estimate an upper limit 
of 2% I-terminated regions within our probing region for 
the data shown here. 

Density-functional theory calculations: Density- 
functional theory (DFT) was performed within the gener¬ 
alized gradient approximation (GGA). We used the DFT 
code, OpenMX [36], based on the linear-combination- 
of-pseudo-atomic-orbitals (LCPAO) method [37] . Spin- 
orbit interaction was included via the norm-conserving, 
fully relativistic j-dependent pseudopotential scheme in 
the non-collinear DFT formalism [38]. We model the Te- 
terminated BiTel surface electronic structure via a super¬ 
cell calculation including a slab of 60 atomic layers and a 
vacuum region of ^20 A thickness. To calculate the spin 
and orbital angular momentum for specific ^-points, we 
used the LCPAO coefficients of local atoms. 
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Fig. SI: Photon energy-dependent ARPES 

Fig. S2: Calculated surface electronic structure of BiTel 

Fig. S3: Coupled Te spin-orbital texture 
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FIG. SI. Photon energy-dependent ARPES. (a) MDC at the Fermi level (Ef d= 15 meV) measured along the Y — M 
direction as a function of photon energy using p-polarised light. This shows non-dispersive features indicative of two-dimensional 
electronic states, (b) This is confirmed by the Fermi momenta extracted by fitting MDCs. The straight lines show the average 
extracted for each band. We thus assign all eight Fermi crossings as arising from two-dimensional subband states formed 
within a near-surface quantum well arising from the downward band bending potential at the Te-terminated surface. This is 
consistent with their spectral weight distribution: the states arising from the second subband (SB2), being less strongly bound 
than SB1 within the almost-triangular quantum well potential at the surface, will have an envelope wavefunction with a larger 
spatial extent perpendicular to the sample surface [39]. This would be expected to give rise to a periodic variation of the 
photoemission spectral weight with photon energy peaked over a narrower range than for SB1 [40], as observed experimentally. 
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FIG. S2. Calculated surface electronic structure of BiTel. (a,d) Atomic weight and (b,e) ratio of out-of-plane to 
average in-plane orbital makeup of the surface electronic structure of BiTel, from Te-terminated supercell calculations within 
density-functional theory including spin-orbit coupling. (c,f) Equivalent orbital-projected electronic structure from calculations 
neglecting spin-orbit coupling. Projections onto the (a-c) first and (d-f) second Te-Bi-I trilayers of the supercell are shown. 
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FIG. S3. Coupled Te spin-orbital texture. Spin textures for Te (a) (b) p y , and (c) p z projections of CESs at 200 meV 

above the Dirac point. The in-plane spin texture is shown by arrows, the out-of-plane by the background colour, both in units 
of h/2. (d) Schematic of the in-plane spin textures coupled to the net in-plane and out-of-plane orbital textures. 






